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ABSTRACT 

The milhmeter-wave rotational emission hnes (4o4 — 3o3 and 5o5 — 4o4) of 
protonated carbon dioxide, HCO^(HOCO"'"), has been detected toward the 
low-mass class protostar IRAS 04368+2557 in LI 527 with the IRAM 30 m 
telescope. This is the first detection of HCO^ except for the Galactic Center 
clouds. The column density of HCO^ averaged over the beam size (29") is 
determined to be 7.6 x lO^'' cm~^, assuming the rotational temperature of 
12.3 K. The fractional abundance of gaseous CO2 relative to H2 is estimated 
from the column density of HCO^ with an aid of a simplified chemical model. 
If the HCO^ emission only comes from the evaporation region of CO2 near 
the protostar (T>50 K), the fractional abundance of CO2 is estimated to be 
higher than 6.6 x 10^^. This is comparable to the elemental abundance of 
carbon in interstellar clouds, and hence, the direct evaporation of CO2 from 
dust grain is unrealistic as a source of gaseous CO2 in L1527. A narrow line 
width of HCO J also supports this. On the other hand, the fractional abundance 
of CO2 is estimated to be 2.9 x 10"'', if the source size is comparable to the 
beam size. These results indicate that gaseous CO2 is abundant even in the 
low-mass star-forming region. Possible production mechanisms of gaseous CO2 
are discussed. 



Subject headings: ISM: molecules, ISM: individual(L1527) 



1. INTRODUCTION 

Carbon dioxide (CO2) is an abundant and important constituent in planetary atmo- 
spheres as well as cometary comae, and is considered to be a key species which links between 
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interstellar chemistry and planetary chemistry. So far solid CO2 was found ubiquitously to- 
ward intermediate- to high mass star-forming regions through observations of its vibrational 
bands with Infrared Space Observatory (ISO) and Spitzer Space Telescope (e.g. de Graauw 
et al. 1996; Grakines et al. 1999; Nummelin et al. 2001; Sonnentrucker et al. 2006). A 
typical abundance of solid CO2 relative to H2 is as high as 10~^ to 10~^ (Gerakines et al. 
1999). Furthermore, Dartoris et al. (2005) discovered huge amount of solid CO2 even toward 
the low-mass class protostar, L723, by observations with Spitzer Space Telescope. Gaseous 
CO2 was also detected in hot regions around several massive young stellar objects, where 
a typical abundance relative to H2 was reported to be an order of 10^^ (e.g. van Dishoeck 
et al. 1996; Boonman et al. 2003b; Sonnentrucker et al. 2006). These observations clearly 
indicate that CO2 is a major chemical component in molecular clouds. 

However, the formation processes of CO2 are not well understood. In particular, a 
contribution of the gas phase production of CO2 is quite uncertain. Although the gas- 
phase chemical models predict the CO2 abundance relative to II2 of about 10"^ (e.g. Lee 
et al. 1997), this has not been confirmed observationally in cold molecular clouds and low- 
mass star-forming regions. A difficulty in measuring the abundances of gaseous CO2 in 
these regions comes from a lack of a permanent dipole moment of CO2. Because of this, 
radio observations of the rotational emission lines are impossible. Infrared observations of 
its vibration-rotation transitions would also be difficult because of insufficient brightness of 
background infrared sources or even of a lack of the sources. An alternative way is to observe 
the rotational emission lines of protonated carbon dioxide, HCO^, in the millimeter- wave 
region, because its abundance is directly related to the abundance of CO2. However, HCO^ 
has never been detected in any cold molecular clouds and any star-forming regions except 
for the Galactic Center clouds like Sgr B2 and Sgr A (Thaddeus et al. 1981; DeFrees et al. 
1982; Minh et al. 1991). It has been thought that HCO^ is produced only under a specific 
condition in the Galactic Center region, where shock chemistry would play an important role 
(e.g. Minh et al. 1988; Charnley et al. 2000). 

An important breakthrough to this situation came from our preliminary line survey 
toward a low-mass protostar IRAS 04368+2557 in L1527 with the Nobeyama 45 m radio 
telescop^ (NRO 45 m), by which peculiar carbon-chain chemistry in this source (Sakai et 
al. 2007, 2008a, 2008b) was being explored. With a very sensitive observation, we recently 
found a weak emission line at the frequency of the 4o4 — 3o3 line of HCO^. IRAS 04368+2557 
is a low-mass protostar in a transient phase from class to class I, which has extensively 
been studied by many researchers (e.g. Ohashi et al. 1997; Hogerheijde et al. 1997, 1998). 
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Hence, this discovery, if confirmed, is very important for understanding of the behavior of 
CO2 during chemical evolution from molecular clouds to protoplanetary disks. Then we car- 
ried out more sensitive observations with the IRAM 30 m telescope (IRAM 30 m) and firmly 
confirmed our detection of HCO^. In this Letter, we report the first detection of HCO^ in 
the star-forming region, L1527, and discuss its astrochemical implication. 



2. OBSERVATIONS 

Observations of the HCO^ lines toward IRAS 04368+2557 in L1527 were carried out 
with NRG 45 m and IRAM 30 m in 2007. We observed the 4o4 - 3o3 and the 5o5 - 4o4 
lines, whose frequencies arc listed in Table 1. The observed position was (q;2000) ^2000) = 
(04^39°^53^89, 26°03'11".0). 

In the observation with NRG 45 m, we used two SIS mixer receivers (S80 and SlOO) si- 
multaneously, whose typical system temperature was about 300 K. The main beam efficiency 
and the beam size of the telescope are 0.51 and 20", respectively, at 85.5 GHz. The telescope 
pointing was checked by observing the nearby SiG maser source (NML Tau) every hour. The 
pointing accuracy was better than 8". The position-switching mode was employed for the 
observations, where the off position was taken at Act = 30', AS — 30'. The backend was 
an acousto-optical radiospectrometer (AOS-W), whose bandwidth is 250 MHz. The velocity 
resolution is 0.87 km s~^ at 85.5 GHz. The intensity scale was cahbrated by the chopper 
wheel method. 

In IRAM 30 m observations, the 3 mm SIS receivers (AlOO and BlOO) were used as 
front ends, whose system noise temperatures ranged from 75 K to 105 K for AlOO and from 
90 K to 130 K for BlOO. The beam size and the main beam efficiency are 29" and 0.78, 

respectively, at 85.5 GHz. The telescope pointing was checked every hour by observing 
nearby continuum sources, and was maintained to be better than 4". The backend was an 
autocorrelator, VESPA. We set the bandwidth and resolution of the individual window to 
be 20 MHz and 20 kHz, respectively. The frequency resolution corresponds to the velocity 
resolution of 0.07 km at 85.5 GHz. The observation was made in the frequency switching 
mode with a frequency offset of 2 MHz. 
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3. RESULTS 

The 4o4 — 3o3 line of HCO^ was first recognized toward L1527 during our preliminary 
line survey with NRO 45 m. Although the line width is slightly broad due to the insufficient 
velocity resolution, the line was detected with the 5.0 a confidence level in integrated intensity 
(Table 1). However, it was still necessary to confirm this detection by observing multiple 
rotational transitions. Then, we made further observations with IRAM 30 m. As a result, 
we successfully detected the 4o4 — Sos and 5o5 — 4o4 lines. Although the signal-to-noise ratio 
of the 5o5 — 4o4 line is not very good, it was detected with 4.7 a confidence level in integrated 
intensity. The line parameters derived from Gaussian fitting are summarized in Table 1, and 
the line profiles are shown in Figure 1. The Vlsr values are consistent with those of other 
molecules in this region (e.g. Sakai et al. 2007, 2008), no significant velocity shift being 
observed. The line widths are slightly narrower than those of the other molecules (~ 0.5 km 
s~^). Because of the narrow line width and sparse spectral density, accidental matching of 
other lines is unhkely (cf. Sakai et al. 2008b). 

The beam averaged column density of HCO^ is determined from the observed intensi- 
ties by the least-squares method. We assume the LTE (local thermodynamic equilibrium) 
condition and take the effect of the optical depth into account. Details of the procedure 
are described elsewhere (Sakai et al. 2008a). The LTE assumption is justified because the 
density of the L1527 core (>10^ cm~^ for ~ 40" scale) (Sakai et al. 2008a) is higher than the 
critical density of the observed HCO^ transitions ((3 — 6) x 10^ cm~^), where the coUisional 
cross sections are taken from Hammami et al. (2007). Since the upper state energies of the 
two observed transitions are not very different from each other, it is difficult to determine 
the rotational temperature of HCO^ signiffcantly. Therefore, we assume the rotational tem- 
perature of 12.3 K, which is the typical value in L1527 (e.g. Sakai et al. 2008a). The column 
density of HCO2" is determined to be 7.6 x VS^^ cm~^, where we employ the dipole moment 
along the a-axis of 2.85 D calculated by the B3LYP/6-311G(d, p) method (Y. Osamura, 
2008, private communication). Even when the assumed excitation temperature is changed 
to 30 K, the derived column density increases only by a factor of 1.8. The optical depth is 
0.006 and 0.005 for the 4o4 — 3o3 and 5o5 — 4o4 lines, respectively, if the beam dilution effect 
is not taken into account. 

The source size of HCOj is unknown, but we could put a constraint on it from our 
observations. The integrated intensity (Tmb) of the 4o4 — 3o3 fine observed with NRO 45 m 
seems to be higher than that observed with IRAM 30 m in spite of a large uncertainty of the 
NRO 45 m data (Table 1). This implies that the HCO^ distribution is at least smaller than 
the beam size of IRAM 30 m (29"). Therefore the peak column density of HCO^ would be 
higher than the beam averaged value obtained above. 
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The present result shows the existence of HCO2 in the low-mass star-forming core. The 
column density of HCO^ in L1527 is less than those toward Sgr B2 and Sgr A (e.g. Minh et 
al. 1988; Deguchi et al. 2006) by almost three orders of magnitude. On the other hand, the 
HCO^ lines were not detected in other regions including star-forming cores such as Orion 
KL and B335 in spite of extensive searches (Minh et al. 1988). One possible explanation 
of the none-detection would be insufficient sensitivity of the previous observations. In fact, 
the column density of HCOj in L1527 is lower than the upper hmits reported by Minh et 
al. (1988) (~ 10^2 cni-2). 



4. DISCUSSION 

4.1. Fractional Abundance of CO2 

In this section, we derive the fractional abundance of the gaseous CO2 with an aid of a 
simplified chemical model. The abundance of HCO^ is related to that of CO2 through the 
following reactions; 

H+ + CO2 ^ HCO+ + H2 : (1) 

HCO+ + CO ^ HCO+ + CO2, (2) 

where H3 is formed by the cosmic-ray ionization of H2 and subsequent reaction of -l- H2, 
and is destructed by 

H+ + CO ^ HCO+ + H2. (3) 

The electron recombination is less important as destruction processes of HCO2" and H^, 
and is ignored. When we assume the steady-state condition, ratio of the number density, 
n(HC0^)/n(C02), can be written as 

n(HCOt) kM^2) si 
n(C02) - hhui^Of /8o^(H2)' 

where fci, A;2, and ^3 represent the rate coefficients for reactions (1), (2), and (3), respectively, 
and C represents the cosmic ray ionization rate, 1.3 x 10~^^ s~^ (e.g. Aikawa et al. 2008). 
Here the unit of n(X) is cm~^, and the fractional abundance of CO relative to H2 is denoted as 
/co- We employ the ki, k2, and ks values of 1.9 x 10"^ cm=^ s"^ (Burt et al. 1970), 7.8 x 10'^° 
cm^ s"^ (UMIST data base), and 1.7 x 10"^ cm^ s"^ (Kim et al. 1975), respectively. The 
numerical factor in equation (4) could be larger or smaller by a factor of 3. 

It should be noted that the n(HC0j)/n(C02) ratio is inversely proportional to the H2 
density, as far as /co is constant. In other words, the fractional abundance of CO2 relative 
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to H2 (/CO2) is simply proportional to n(HC02 ); 



/C02 - 



5.4 X 10^ X /^o X ^(HCO^). 



(5) 



If CO2 exists only within the sphere with a radius of R around the protostar, with the 
constant /co2) ^^(HCOg ) is also constant and is related to the beam averaged column density 



assuming the Gaussian-shaped telescope beam with an e-fold beam radius of Ri,. Note 
that the right-hand side of equation (6) approaches to 2i?n(HC0^) for R ^ R^, and to 
n(HC02 )(47ri?^/3)/(7rii!^) for it! ^ R^. With equation (6), we can derive /co2 directly from 

the observed A^ave(HCOj) without any assumptions on the H2 density and the H2 column 
density We calculate /co2 as a function of i?, as shown in Figure 2, where we employ the 
/co value of 3.9 x 10~^ reported by J0rgensen et al. (2002). The half-power beam radius 
{RbVlri2) for IRAM 30 m is 2000 AU, where the distance to L1527 is assumed to be 140 pc. 

According to Marct ct al. (2004), the size of the hot region where CO2 can be evaporated 
(>50 K) is estimated to be smaller than 140 AU for L1527. If gaseous CO2 traced by 
our HCO^ observation are confined to this small radius, fco-z should be evaluated to be 
higher than 6.6 x 10"'* by using Figure 2. In the ice evaporation rcgion(>100 K), H2O and 
NH3 would become very abundant and contribute to further destruction of HCO2 in the 
innermost region. This effect can be involved as an effectively larger /co value in equations 
(4) and (5). Then /coa is comparable to the elemental abundance of carbon in interstellar 
clouds (7 X 10~^). On the other hand, Furlan et al. (2007) recently reported the mass 
fractional abundance of solid CO2 relative to the gas in L1527 to be 1 X 10"^ with the 
Spitzer observation, which corresponds to the fractional abundance of solid CO2 relative to 
H2 of 5 X 10^®. This is much lower than the /co2 value derived from our observation. Hence, 
the direct evaporation of CO2 is very unrealistic as a major source of gaseous CO2 traced 
by HCO^ in L1527. This would further be supported by the narrow line width (0.57 km 
s~^) of HCO^ in comparison with the line width of the very high excitation line of HC5N 
(J = 32 - 31, 85.2 GHz, E^^Gl K) observed with IRAM 30 m in L1527 (Sakai et al. in 
preparation) . 

If we assume /co2 of ~ 5 x 10^^, which is close to the value reported in Orion-KL IRc2 
by the infrared observation (Boonman et al. 2003a), R is derived to be 1600 AU. On the 
other hand, /co2 is still as high as 2.9 x 10"^, even if HCO^ is extended over the beam size 
of IRAM 30 m. Therefore, substantial amount of CO2 exists even in the low temperature 
region where thermal evaporation of CO2 from grain mantles is very inefficient. It should be 
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noted that the interaction between outflows and ambient gas would not be responsible for 
the CO2 production in L1527, because the line width of HCO2 is narrow and no velocity 
shift is observed. 



4.2. Origin of Gaseous CO2 

As described above, most of gaseous CO2 in L1527 does not originate from direct evap- 
oration of solid CO2, but from gas- phase reactions. Since the gas phase production of CO2 
is not well understood observationally (e.g. van Dishoeck et al. 2004), we discuss a few 
possibilities. First we compare our result with the chemical model calculation by Lee et 
al. (1997). According to their 'new standard' model, the gas phase abundance of CO2 is 
expected to be 2.8 x 10^^ and 1.4 x 10^^ at 10^'^ yr and the steady state, respectively, for 
the H2 density of 10^ cm~^. This is almost consistent with our result in L1527, if the CO2 
distribution is extended over the beam size of IRAM 30 m. 

Secondly, CO2 might be produced in a hot region through a reaction of OH and CO. 
Since this reaction has activation barrier of 176 K (UMIST data base; See also Talbi and 
Herbst 2002 for the potential energy surface.), it is efficient in a small region near the 
protostar, like a hot core or a hot corino. In such a hot region, the abundance of H2CO 
is enhanced due to mantle evaporation, which would supply HCO through the electron 
recombination reaction of H2COH+ formed from H2CO and H3 . Then, CO2 will be produced 
through the HCO -|- O reaction (e.g. Snyder et al. 1985). However, these high temperature 
mechanisms occur in a much smaller region than the CO2 evaporation region, and hence, 
they can not be major contributors. 

Since L1527 shows warm carbon-chain chemistry (WCCC) (Sakai et al. 2008a), we 
finally consider a possible CO2 formation related to the WCCC. In the WCCC, CH4 evapo- 
rated from grain mantles drives the carbon-chain chemistry. At the same time, HCO would 
be produced by various reactions, which successively forms CO2 by the reaction with the 
atomic oxygen. The warm region, where CH4 can be evaporated, is substantially larger than 
the region where CO2 can be evaporated, and hence, CO2 would be distributed over the size 
of WCCC (Sakai et al. 2008a). Rather compact distribution of HCO^ inferred from the 
beam dilution effect, if correct, favors this scenario rather than the production of CO2 in a 
cold cloud. In this relation, the HCO (Iqi — Oqo, 86.7 GHz) fine is detected toward L1527 
in our preliminary line survey, which may also support the above idea. Furthermore, an 
enhancement of gaseous CO2 in the warm region, where CH4 is evaporated, is predicted in a 
chemical model of a dynamically evolving cloud by Aikawa et al. (2008, Figs. 6 and 10). In 
order to establish the formation pathway of CO2, distribution of HCO2 in L1527 is a key. A 
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high spatial resolution observation with an interferometer would be necessary. By comparing 
the result with the detailed chemical model calculations, we will be able to investigate the 
contributions of cold gas phase chemistry and WCCC to the formation of CO2. Sensitive 
observations of HCOj toward various clouds would also be needed to understand the origin 
of gaseous CO2. 

We are grateful to Eric Herbst for his invaluable discussion. We also thank to Yoshihiro 
Osamura for communicating the results of his ab intio calculations. We thank to the staff 
of the IRAM 30 m telescope and the NRO 45 m telescope for their excellent support. This 
study is supported by Grant-in- Aids from Ministry of Education, Culture, Sports, Science, 
and Technologies (14204013, 15071201, and 19-6825). 
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Table 1: Observed Line Parameters for HCOg 



Transition 


Frequency 




T/a 


rms*^ 


/ TMBdv{3a) 




[GHz] 


[K] [km s-i] 


[km s-1] 


[mK] 


[K km s-i] 


IRAM 30 m 


4oi — 3 03 


85.5315123 


0.058(8) 0.32(5) 


5.88(2) 


4.5 


0.020(4) 




106.9135639 


0.044(6) 0.42(6) 


5.88(3) 


10.1 


0.020(13) 


NRO 45 m'^ 


4o4 — 3o3 


85.5315123 


0.039(11) 1.0(3) 


6.2(2) 


8.0 


0.040(24) 



Note. — The numbers in parentheses represent the errors in units of the last significant digits. 

^ Obtained by the Gaussian fit. 

^ The rms noise averaged over the linewidth. 

Five successive channels of VESPA are bound to improve the S/N ratio. 

The line width becomes broad due to the insufficient velocity resolution of AOS-W. 
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Fig. 1. — Spectral line profiles of HCO^ (4o4 — 3o3 and 5o5 — 4o4) observed toward L1527. 
Negative features at -1 and 13 km s~^ in the 4o4 — 3o3 spectrum are frequency-switching 
artifacts. 



- 12 - 




Fig. 2. — Fractional abundance of CO2 relative to H2 as a function of the assumed radius of 
the emitting region of the HCO^ line. 



